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Abstract 

We construct quartic quasitopological gravity, a theory of gravity containing terms quartic in 
the curvature that yields second order differential equations in the spherically symmetric case. Up 
to a term proportional to the quartic term in Lovelock gravity we find a unique solution for this 
quartic case, valid in any dimensionality larger than 4 except 8. This case is the highest degree of 
curvature coupling for which explicit black hole solutions can be constructed, and we obtain and 
analyze the various black hole solutions that emerge from the field equations in (n + 1) dimensions. 
We discuss the thermodynamics of these black holes and compute their entropy as a function of the 
horizon radius. We then make some general remarks about K-ih order quasitopological gravity, 
and point out that the basic structure of the solutions will be the same in any dimensionality for 
general K apart from particular cases. 
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I. INTRODUCTION 



The gauge/gravity idea is that gravitational dynamics in a given dimensionality can be 
mapped onto some other (nongravitational) field theory of a lower dimensionality The du- 
ality between central charges and couplings on the nongravitational side and the parameters 
on the gravitational side has been explored primarily through the trace anomaly [l]. How- 
ever Einstein gravity does not have enough free parameters to account for the ratios between 
central charges and therefore is only dual to those conformal field theories for which all the 
central charges are equal. To broaden the universality class of dual field theories which 



one can study with holography, one must extend to more general theories of gravity which 
contain more free parameters such as Lovelock theory Q or quasitopological gravity {3, ^]. 
These additional central charges have recently been investigated holographically {5, If]]. 

Another point which is interesting in gauge/gravity duality is that the dual CFT should 
respect causality. This creates a constraint on the coupling constants of the gravity the- 
ory. In this analysis, one considers graviton fluctuations that probe the bulk geometry in 
the presence of a black hole. In general the dual CFT plasma may support superluminal 
signals, and so the gravitational couplings must be constrained so as to avoid the appear- 



ance of such superluminal modes. For Lovelock gravity, whi 
match those arising from requiring positive energy fluxes 



e causality constraints precisely 
8], it has been shown that this 



matching does not appear in general, specifically, for cases where the gravitational equa- 
tions of motion are not second order [9(. However for cubic quasitoplogical gravity there are 
three constraints that arise from requiring positive energy fluxes, which determine the three 
coupling constants. No evidence for causality violation was found once the curvature-cubed 
coupling was chosen consistent with these constraints j^]. 

Motivated by the success of holographic studies of second Ml, and third-order Lovelock 

n n n n n 

gravity |8j, |10] and curvature-cubed, or cubic quasitopological gravity |5|, |6|, , we consider 
here adding a quartic curvature term with a new coupling constant on the gravity side, 
affording exploration of a larger space of field theories. In quasitopological gravity, the 
linearized equations in a black hole background are fourth-order in derivatives and so one 
does not expect causality constraints to match those arising from requiring positive energy 
fluxes. In view of the results for the cubic case j^], the simplest nontrivial case to consider 
is the quartic case; with this new coupling constant, we have four coupling constants and 
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therefore the constraints arising from causality may not match the three constraints arising 
from requiring positive energy fluxes. The first step in such an investigation is to construct 
the quartic theory and analyze its basic properties. We shall consider the more detailed 
considerations of positivity of energy and causality in future work. 

Progress with Lovelock gravity and cubic quasitopological gravity relies on the fact that 
even though this is a higher curvature theory of gravity, the holographic calculations in this 
model are still under control, at least in spherically symmetric settings. This control in turn 
is based on the two facts that the equations of motion are only second-order in derivatives 
(again, for spherical symmetry) and that exact black hole solutions have been constructed. 
Hence we want to introduce a quartic curvature topological gravity for which the equation 
of motions are still second-order and exact black hole solutions can be constructed. In this 
context, fourth-order Lovelock and quasitoplogical gravity is the largest order for which the 
field equations can be solved analytically. Even in Lovelock gravity, this largest analytic 
solution has not yet been considered. Here we carry out the first steps along these lines, 
studying the exact spherically symmetric solutions and their properties. 

To obtain the quartic natural generalization would be to add interaction terms 

quartic in curvature via fourth-order Lovelock gravity. However, because of the topological 
origin of the Lovelock terms the quartic interaction term of Lovelock gravity only contributes 
to the equations of motion when the bulk dimension is nine or greater. In the context of 
gauge/gravity duality, this means that such a term will be effective in expanding the class of 
dual field theory in eight or more dimensions. Our key result in this paper is to construct a 
new gravitational action with quartic curvature interactions (quartic quasitopological grav- 
ity) valid in lower dimensions, thereby providing a useful toy model to study a broader class 
of four (and higher) dimensional CFTs, involving four independent parameters. 

Here we explicitly construct quartic quasitopological gravity in any dimensionality except 
8, beginning with the five dimensional case. Although an action quartic in curvature terms 
has been previously constructed [3] (and from which was proved a generalized Birkhoff theo- 
rem, namely that constant spherical/planar/hyperbolic transverse curvature implies staticity 



12(|), the field equations in the spherically symmetric case vanish in less than seven space- 



time dimensions. In contrast to this, the quartic topological action we construct yields 
nontrivial second-order field equations in all spacetime dimensionalities but 8. Indeed, our 
quartic curvature action differs from that of ref. in terms of its various coefficients, and 
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insofar as it yields nontrivial field equations in 5 dimensions and higher. 

We also present and discuss exact black hole solutions of this new theory for various 
asymptotic boundary conditions. These solutions share a number of features in common 
with solutions from higher-order Lovelock theories in greater dimensions. For example, in 
the spherically symmetric cases we consider, the field equations for the metric function 
in our quartic theory in 5 dimensions are formally the same as for fourth-order Lovelock 
theory in 9 dimensions, differing only by the power of r present in the resultant quartic 
equation. We furthermore consider the thermodynamic behavior of these objects for general 
dimensionality. We leave a detailed study of the properties of the dual class of field theories 
for future investigation. 



II. QUARTIC TOPOLOGICAL ACTION IN FIVE DIMENSIONS 

Motivated by considerations of the AdS/CFT correspondence, we want to consider a 
curvature-quartic theory of gravity in five dimensions. We are interested in a gravity theory 
which produces second-order equation of motion and can have exact solutions. A natural 
candidate that has these properties is the fourth-order Lovelock gravity with action 

Jg = / dn+lx ^i~ 2A + £ i + ^2 + a 3 C 3 + « 4 A], (1) 

where A = —n(n — 1)/ 2Z 2 is the cosmological constant Jbr AdS solutions, and the ctj's are 
Lovelock coefficients with dimensions (length) 2l ~ 2 and 



[" . — J_AW2 'M2i D ^1^2 ...J) ^2i-1^2i (n\ 

1 2 i 1 l ' 2 "' L ' 2i WW V2i-lM2i " v / 



A key property of this action is that the term proportional to contributes to the equations 
of motion in dimensions with n > 2k. Hence the above action with interaction terms 
quartic in the curvature tensor contribute to the equations of motion only in nine and 
higher dimensions and hence will not contribute in the desired five dimensions. 

While Lovelock's Lagrangian yields second-order equations of motion for an arbitrary 
spacetime, we limit ourselves to the case of spherically symmetric spacetimes. The metric 
of 5-dimensional spherically symmetric spacetime may be written as 

dv 2 

ds * = -N^ r )f(r)dt 2 + + r 2 d^ 3 , (3) 



4 



where oE| 3 represents the metric of a 3- dimensional hypersurface with constant curvature 
Qk and volume V 3 . The first three terms in the action ([1]) contribute to the field equation in 
five dimensions, while the third and fourth-order Lovelock terms do not. 

Our aim is to include in the action terms quartic in the curvature that contribute to the 
field equations in five dimensions and yield second-order equations of motion for spherically 
symmetric spacetimes. We find that this action may be written as 

= T7T- [ d n+1 x^[-2A + A + /i 2 £ 2 + fx 3 X 3 + 1*4X4], (4) 

107T J 

where C\ — R is just the Einstein-Hilbert Lagrangian, £ 2 = RabcdR abcd — 4:R ab R ab + R 2 
is the second-order Lovelock (Gauss- Bonnet) Lagrangian, and X 3 is the curvature-cubed 
Lagrangian 



V _ D nbedf p a c , * / J l JU ° J p jjabcd p 

X 3 - RabcdR R e f + J7T 7TT- oV o RabcdR H 



1 f3(3n-5) 
(2n - l){n - 3) 
-3(n - l)R abcd R abc e R de + 3(n + l)R abcd R ac R bd 

+ 6(n - l)R ab R bc Rc a - 3{3n ~ 1] R a bR ab R + ^±lii? 3 ) (5) 



obtained previously {4]. 

In eight dimensions, there are 26 distinct scalar functions that are quartic in the curvature 
tensor [13j. However, one may construct the fourth-order Lagrangian of Lovelock gravity by 
combining the following 25 terms 

p4 p2 p pa6 p2 p pabcd 

ti , ti tiahti , ti Habcd-ti , 

ppapbpc pp «p fiD oi p p a p be p de 

-K-ftft rt c rt a , rtrt c rt^ .K ao , rtrtb rtrf e rt ac , 

p p ab p cd p e/ p p ab p cd p e/ P a P^P c P^ 

-K-tt c d -Ke/ -Kab , -K-Kce ^-a/ -Kbd , -Kfe -Kd -"-c > 

papbpcpd papbpcpde p a p b p cd p e/ 

-ttfc rt c rtrf r£ a , -K itrf it e r£ ac , rt;, r£ a r£ e j r£ c rf , 

p a p b p cd p e f p a p 6 p cd p ef p a p b p cd p e f 

-Kfe -K c /t e / ti ad , tic tid ti e f ti ah , ti c ti e ti a f ti hd , 

p a p b p cd p e/ p a p fee p de p /g p a p 6c p de p /g 

-n-c -Ke -n-6/ -Kad , -Kb ti a d ttfg ti ce , it;, tide tLfg ti ac , 

pap 6c p de p /a / pafced p \2 pa6cd p ejyfgh p 

-K& -Kd/ -flag -K C e , -KabcdJ , -f<abc d-tifghe, 

ryabed p pefg/ip p ab p cd p ef p g/i p ab p cd p ef p g/i 

-K ti e f cd ti J » tiabgh, tied tt eg ti ah J ti h f , ti ce ti nn ti hh ti,]f , 

p ab p cd p ef p g/i 
-Kce "^9 K dh ti h f . 

in a particular way 
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Since the X3 Lagrangian contains no derivatives of the curvature tensor, we shall construct 
the X4 term in the action using only the above 25 terms. For the metric ([3]) the function 
N(r) performs the role of the lapse function, making it possible to write the action as a 
functional of f(r) and its derivatives, with N(r) appearing linearly in the action. 

Since the Riemann tensor has at most 2 derivatives of the metric functions we find that 
there are at most 8 derivatives in any term for the quartic-curvature action. We require 
all terms in the Lagrangian to vanish that have more than two derivatives. For the metric 
ansatz ([3]) not all 25 terms above are needed to ensure the resultant equations of motion are 
second-order differential equations. Remarkably we find that we can choose 

X 4 =ciR a t )Cl iR cde f R hg e fRh g ab + C2R a bcdR abcd RefR e ^ + c 3 RR a f,R ac R c b + C4(R a t, c dR abcd ) 2 

I ^ D P ac D T)db I „ ~D ~D P ac T>db i « D D ac D& e T}d , p pace/ pb pa! 
+C5KabK rC cc irC + CeKK a bcd^i tX + Cjriabcdtt ti ti e + C^tiabcdti It e K j 

+CgR a bcdR aC Ref R^ + Cio-R 4 + CuR? ' R a b c dR abcd + C\ 2 R 2 RabR^ 

+ci3R a bcdR R e f c g R d9 + ci A R abcd R aec f R geh fR gbhd , (6) 

without loss of generality. We must then choose the coefficients q to yield only a second- 
order contribution to the field equations. We find that 

ci = -1404, c 2 = 1848, c 3 = -25536, c 4 = -7422, c 5 = 24672, 
c 6 = -5472, c 7 = 77184, c 8 = -85824, c 9 = -41472, c w = -690, 
en = 1788, C12 = 6936, c 13 = 7296, c u = 42480 (7) 

is the unique solution up to a term proportional to the quartic Lovelock Lagrangian. 
Defining the dimensionless parameters fa .... fa to he 

I 2 2 4 21024 

fa = — gA, fa = -/i 2 , fa = fa = —j^Vi (8) 

and integrating by parts, we find that the action (TjJ per unit volume V3 reduces to the rather 
simple form 

3 



dtdrN(r) {r 4 (fa + ip + faip 2 + faip 3 + fa^) }' , (9) 



" 16tt/ 2 „ 

where prime denotes the derivative with respect to r and ip = Pr~ 2 (k — f). 
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III. GENERALIZATION TO n + 1 DIMENSIONS 



In this section we consider the action (j3J) in n + 1 dimensions for the spherical metric 

dv 2 

ds 2 = -N 2 (r)f(r)dt 2 + r 2 dY? Kn ^ (10) 

fir) 

where dT? kn _ 1 represents the metric of an (n — l)-dimensional hypersurface with constant 
curvature (n — l)(n — 2)k and volume V n -\. Using the same procedure as in the preceding 
section for five dimensions, we can obtain the coefficients Cj's in Eq. (jH]). The results are 
somewhat cumbersome so we list them in the Appendix. 

As before, we find that after integrating by parts and defining the dimensionless param- 
eters /to, £12, A3 and /t4 to be 

2l 2 A (n-2)(n-3) A (n - 2){n - 5)(3n 2 - 9n + 4) 

^ ~ ~n(n^T) ^ " P M2) ^ " 8(2n - 1)1* ^ 

A ^ n (n - 1) (n - 2) 2 (n - 3) (n - 7) (n 5 - 15 n 4 + 72 n 3 - 156 n 2 + 150 n - 42) 

/i4 — 7g ^4, 

the action per unit volume reduces to 

Ig = / dtdrN(r) [r n (fi + if; + fi 2 if; 2 + fa^ + M*)]' , (11) 

where again if) = l 2 r~ 2 (k — /). Note that in the absence of a cosmological constant /to = 0, 
while in the presence of a positive/negative cosmological constant we take /to = ±1. 

We pause to comment that /t 4 is zero in 8 dimensions, suggesting that X± yields another 
topological invariant in eight dimensions besides the 8-dimensional Euler density (given by 
L4 in Eq. fl2])). However it is straightforward to show that A4 has 8-th order derivative 
terms for nontrivial 8-dimensional geometries and therefore is not a topological invariant. 
Hence we refer to this theory of gravity as quartic quasitopological gravity. Note that our 
construction does not yield a nontrivial quartic interaction term in n + 1 < 4 as well. 

Varying the action (jUj) with respect to N(r), we obtain 

[r*(/2a + V + £2^ + M 3 + M*)]' = (12) 

for the equations of motion. Formally this equation is the same as that obtained from 9- 
dimensional fourth-order Lovelock gravity in the spherically symmetric case. However the 
power of r differs in (fl2|) from this case, since n > 5 can have any integer value except 8. 
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The black hole solutions to this equation will consequently have analogous properties. They 
will be asymptotically flat, AdS, or dS depending on the choice of parameters (as we will 
discuss below) and they will have a scalar curvature singularity at r = cloaked by an event 
horizon. 

The solutions (1121) are the real roots of the following quartic equation 

^ + + + L$ + J_« = o, (13) 

A^4 A^4 A^4 A^4 



where 



TYl 

« = &-—, (14) 



and m is an integration constant which is related to the mass of the spacetime. 

The geometrical mass of black hole solutions is 

ft 2 I 4 . I 6 I 8 \ 

m = I + k— + fak 2 — + /U3^ 3 — + frik 4 — ) r+ (15) 
\ r + r + r + r +/ 

in terms of the horizon radius r+. Before considering the properties of particular solutions, 

we compute the Hawking temperature 

_ 1 n/z r+ + (n - 2) A;/ 2 r^ + ( n - 4) k 2 fi 2 l 4 r\ + (n - 6) k 3 fi 3 l 6 r 2 + + (n - 8) A; 4 /^/ 8 

~4 ( + 2kfi 2 l 2 r 4 + + 3k 2 fi 3 l 4 r 2 + + 4/} 4 A; 3 Z 6 ) vr / 2 r + ^ 

for the general black hole solution given by Eq. (fl9|) . Clearly, T always positive for k = 0, 

and therefore there is no extreme black hole. However, for k — ±1, extremal black hole 

solutions exist with horizon radius r ext , where r ext is the largest real root of 

nfar 8 ext + (n - 2) A;/ 2 r e 6 xt + (n - 4) fc 2 /i 2 / 4 r 4 xt + (n - 6) Wr^t + (n - 8) A: 4 /^/ 8 . (17) 

Equation (TTTj) has at least one real solution in the absence of a cosmological constant (/to = 
0). Hence there exist black holes with inner and outer horizons, extreme black holes or 
naked singularities, depending on the choice of parameters. 

However for nonzero cosmological constant, extreme black holes appear as solutions pro- 
vided 

B 2 

A = A 3 + — > 0, 

where 



A = 3 (n — 2) (n — 6) jl 3 — Yin (n — 8) /to/t4 — (^ — 4) /t 2 , 
5 = -9 (ra - 4) fa [8n (n - 8) fafa + (n - 2) (n - 6) /t 3 ] 

+27(n - 2) 2 (n - 8) /2 4 + 27n (n - 6) 2 /to/tg + 2 (n - 4) 3 fi\. 

8 



The mass of the extreme black hole may be obtained by using Eq. (fl5l) and computing 
m ex t — iTiifext)- Then, our solution corresponds to a black hole with inner and outer horizons 
provided m > m ext , an extreme black hole if m = m ext , and a naked singularity for m < m ext 



IV. SPECIAL SOLUTIONS 



We consider in this section special solutions of Eq. f fl3|) . Eliminating the cubic term in 
Eq. ( fl3|) by use of the transformation 



yields 



where 



* = (18) 



X 4 - a X 2 + Px ~ 1 = 0, (19) 



S/t4 jj,4 ' 8/i| 2jl1 /t 4 ' 

16(4 /t 4 /X4r n ' 



The most general solution of ffl9|) will yield the most general metric solution for quartic 
quasitopological gravity with constant curvature horizons. However special solutions will 
emerge for particular choices of the coefficients; we first examine these. 

A. a = 0, P = 0: 

In this case (13 and /14 are: 

and the asymptotically AdS solution is 



3 r 2 J / 8// 2 



772 



1/4' 



^w^+^fH 1 3 r-^JJ >■ (21) 

Requiring nonsingular real solutions implies that < fi 2 < 3£i /8. 

For /to = 1, the minus branch corresponds to an asymptotically AdS black hole solution 
with two horizons provided m > m ext , an extreme black hole if m = m ext , and a naked 
singularity for m < m ext (see Fig. [1]). The plus branch always yields a naked singularity 



FIG. 1: The asymptotically anti de Sitter case (A): f(r) vs r for k = 1, n = 4, fa = 1, fa = -2, 
1=1 and m < m ex t, m = m ex t and m > Tn ex t from up to down, respectively. 



for k = 0,1. However for k 
event horizon is located at 



-1 it corresponds to a black hole with a single horizon. The 



2/i 3 

«) ' Yi 



-X 



+ 



2 x n " 4 -ml- n x A 



(22) 



where x = r+/l. In this case the mass parameter can even be negative above a certain lower 



bound 



15) . A similar situation holds for fi = —1. 



Although such solutions do not have a smooth general relativistic limit as fa 0, it 
is possible that phase transitions to this branch from the minus branch could occur. This 

n 

phenomenon has been demonstrated to take place in Gauss-Bonnet gravity [16(. Despite 
both branches having positive mass 1^, the plus branch is perturbatively unstable. Quan- 
tum transitions can occur between the two vacua, and neither the empty Einstein vacuum, 
nor the empty Gauss- Bonnet vacuum provide a good description of the stable quantum 



16) . Whether or not a 



vacuum, since each becomes populated with bubbles of the other 
similar phenomenon takes place in quasitopological gravity remains an interesting topic for 
future investigation. With this in mind, we will henceforth consider only the minus branch 
of the solutions. 
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For /to = — 1 and k = 1 asymptotically de Sitter solutions are present for the minus 
branch. These correspond to black holes with two horizons, an extremal black hole with 
one horizon, or a naked singularity, depending on the relative size of m. Asymptotically 
dS black holes exist provided m ext < m < m crit , where m ext and m cr ; t are the values of 
the mass parameter for the smaller and larger root of T = respectively. We illustrate in 
Fig. [21 the behaviour of the metric function f(r) for the various cases. In the absence of a 




FIG. 2: The asymptotically de Sitter case (A): f(r) vs r for k = 1, n = 5, /to = — 1, fi2 = -4, 
1=1 and m < m ex t, rn = m ex t, m cxt < m < m cr i t ,m = m cr it, and m > m cr ; t from up to down, 
respectively. 



cosmological constant (/to = 0), the case k = 1 yields an asymptotically flat black hole with 
metric function 

8/t 2 m \ 1/4 



f(r) 




1 



(23) 



2/t 2 P 

This solution corresponds to a black hole with two horizons provided m > m ext , an extreme 
black hole if m = m cxt , and a naked singularity for m < m ex t (see Fig. [3]). 
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FIG. 3: The asymptotically flat case (A) :f(r) vs r for k = 1, n = 4, /to = 0, /t2 = -5, 1=1 and 
m < m ext , to = m cxt and m > m cxt from up to down, respectively. 

B. /3 = 0: 



Another special solution of Eq. ffl~9l) corresponds to the case of f3 = 0, for which Eq. ( flUj) 
is quadratic in x 2 an d 

fl'2 



4/t 3 /t 4 



The metric function /(r) can be written as 



/ 



/(r) 



r 

7^ 



^3 



± 



V 



/':; 



16/t| 



4/t 4 ^ I6//3//4 



+ 



/i 3 /i, 4 



+ 



/i 4 r r 



(24) 



/ 



Since we are interested in black hole solutions, we choose the minus branch of /(r) for 
= 0, 1. For /t = 1, the minus branch of this solution corresponds to an asymptotically 
AdS black hole with two horizons, an extreme black hole or a naked singularity provided 
m > m ex t, m = m ext , and m < m ex t, respectively. In Fig. H]we illustrate the various cases. 
For /t = — 1 and k — 1, the solution corresponds to an asymptotically de Sitter black hole 
with two horizons if m ex t < m < m cr ; t , an extremal black hole with one horizon if m = m ex t, 
or a naked singularity otherwise. We illustrate the different possibilities in Fig. [5j 
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FIG. 4: The asymptotically anti de Sitter case (B): f(r) vs r for k = 1, n = 4, /jlq = I, fa = .2, fa = 
.015, 1=1 and m < m ex t, m = m ext and m > m ex t from up to down, respectively. 



For zero cosmological constant and k — 1, the metric function for the asymptotically flat 
black hole solution is shown in Fig. [6] for various mass parameters. For k = — 1 we also find 
black hole solutions with a single horizon for the plus branch; we have not illustrated the 
metric function here. We note, as per our earlier discussion, that this plus branch has no 
smooth Einsteinian limit. 

V. GENERAL SOLUTIONS 

We consider first asymptotically (A)dS solutions, for which we have two real solutions 
provided A at infinity (n — > /to) is positive where 



A 



27 + T 



and 



C 



3/t3 — fi\ 4k 



(25) 



3/^4 /i4 



D 



27 fil 3 Ul 




(26) 
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FIG. 5: The asymptotically de Sitter case (B): /(r) vs r for k = 1, n = 5, /to = —1, A2 = -4, 
//3 = .08, 1=1 and m < m ext , m = m ext , m ext < m < m cr ;t ,?n = m cr i t , and m > m. cr it from up to 
down, respectively. 



The real solutions of Eq. f[T3|) are 

f(r) = 

where 

i2 = 
E = 



r 



J4_ 


2/t 2 




3/i 4 




2/i 2 


4/il 


/t 4 



, t 4 + -i2±^ 
/ 2 V4/k 2 2 



1/3 + (I- 



(27) 



1/3 



1 

4R 



4//2A3 8 



1/2 



1/2 



(28) 
(29) 



describing the two physical branches of the solution. 

Again we are interested in black hole solutions that have a smooth Einsteinian limit. 
Therefore we choose the minus branch of f(r) 

.2 / ,._ l j 



(30) 



Figure [7] shows the metric function f(r) for different values of mass parameters with jj, = +1. 
For k = —1 the solution yields a black hole with one horizon. 



14 







2 



4 



6 



8 



10 



-0.2 



r 



-0.4 



FIG. 6: The asymptotically flat case (B): f(r) vs r for k = 1, n = 4, /to = 0,jj-2 = -2, £13 = .015, 
1=1 and m < m ext , m = m ex t and m > m ext from up to down, respectively. 

For /to = — 1 and fc = 1, the solution is that of an asymptotically de Sitter black hole 
with two horizons if m ext < m < m CT[t , an extremal black hole with one horizon if m = m ext , 
or a naked singularity otherwise (see Fig. [8]). 

Second, we consider asymptotically flat solutions. These are present only if /t = (the 
cosmological constant vanishes), implying that Eq. ( jl~3l) reduces to 



in the large-r limit. We see from Eq. ( I3T1) that one can have an asymptotically flat solution 
with = 0. This asymptotically flat solution can be written down by substituting /t = 
in Eq. (!30|) . It corresponds to a black hole with inner and outer horizons provided m > m exk , 
an extreme black hole if m = m cxt , and a naked singularity for m < m ext where m ext can be 
calculated numerically. The metric functions for these black holes are shown in Fig. [91 




(31) 
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FIG. 7: The general asymptotically anti de Sitter case: f(r) vs r for k = 1, n = 4, fa = 1, fa = -2, 
//3 = .l,//4 = .06, 1=1 and m < m ex t> m = ^ext and m > m ext from up to down, respectively . 



VI. ENTROPY DENSITY 

The entropy of the black hole solutions can be calculated through the use of the formula 

Q 

S = -27T (h d n 1 x^/§ — e ab i cd , 

J Otiabcd 

where £ is the Lagrangian, g is the determinant of the induced metric on the horizon and 
iat is the binormal to the horizon. For the static black holes considered here 

v- J^-~ ~ 

ar> ^ab^-cd 
O-tlabcd 

is constant on the horizon and so the entropy density is s = S/V n -\ = — 2-nr r ^T 1 Y . For 
Einstein Lagrangian, Y\ = — l/(87r) and the resulting entropy density is the expected 
Bekenstein-Hawking entropy s = l/4r" _1 [3]. Applying this formalism to Gauss-Bonnet 
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FIG. 8: The general asymptotically de Sitter case: f(r) vs r for k = 1, n = 5, /Iq = — 1, fi2 = -4, 
/f 3 = .l,/t4 = .002, 1=1 and m < m ext , m = m ext , m ext < m < m crit m = m crit , and m > m crit 
from up to down, respectively. 



and cubic terms, one obtains 



4] 



S2 



2(n-2)(n-3) 
(n- l)/ 2 



(i? - 2 (ij* + BT r ) + 2R\) 



2(n - 3)r 2 



kfx 2 r+ 



n-l 



S 3 



4/x 3 Z 4 r 



A n-l 



(n - 2)(n - 3)(n - 5)(3n 2 - 9n + 4) 
3(n - 1) 



3(n-3)(2n-l) 



(32) 

/ Titm prn r)tm pr n\ 

\ -ft tn-ft rra rra-ft mt J 



jv ( pir p m ptr p m , J_ / p pmnpr , p jjmnpt\ \ 

g I tm-L^r -ft rm-L^t ~T~ ^ ^-fi'trmpr-'l' ~T~ -fLmript -ft ) J 

9(n- 1) 



3(3??- ^) / op ptr ■ p p/m,/«/ 1 

•— — | Z_ft_ft tr + -K mnpq n | + 



(n-3) V 
3(n + l) 



1 

2" 



2(n - 3) 
3(3n - 1) 



2(n-3) 

^-ft t-ft r — IX T 1X t T~ -ft mrnJx + -ft m f n -ft ) 



1^-rt -ft rm + -ft -fttm J 



4(ra - 3) 



(RmnR mn + -R (-R ? r + ^t)) + ~tR^ 



16(n-3) 



3(w- l)/ 4 
4(n — 5)r 4 



7c 2 /f 3 r 



n-l 

+ ' 



(33) 
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FIG. 9: The general asymptotically flat case: f(r) vs r for k = 1, n = 4, /to = 0,/i2 = -2, 
= .l,//4 = .06, 1=1 and m < m ex tj 771 = m ext and W > m ex t from up to down, respectively. 



respectively. We can use the same formalism to obtain the entropy of the quartic term ([6]). 
It is a matter of calculation to show that Y4 reduces to 



16vr ( n - 7)n{n - l)(n - 2) 2 (n - 3)(n 5 - 15 n 4 + 72 n 3 - 156 n 2 + 150 n - 42) 

16CiR rte f R h9 e fRhgrt + 2C2{^R e fR e R tr tr + RabcdR abCd (R* r + + C 3 (2R a bR aC R c b 

+3R(R ar R ar + R at R a t)) + lQc4R a bcdR abcd R tT tr + ^c^{R a b{R a r R br + R a tR bt )) 

+2CQ(R(R t tR 7 r — R t rR r t) + RR aC (Ra cr + Ra ct) + RabcdR^R^" 1 ) 

+C7^rt j-it rite + ri t ri K re — ti r ti ti te — ti t ti ri re + ti ti e [K brd + ti btd) 

1 op pac pdr . op pac pc(t\ . / p re/ p pt i p te/ p pr O D' re l D E> 
+^-Karcd-K -K + ^-tiatcd-ti ti ) + C 8 {ti r J ti te ti f + ti t ti re ti f — IK J ti re ti t f 

1 o p / p6r pcit p6i pdr \ , o pd / p parce , p paice \ \ 
+ZK r btd{ti ti — ti ti ) + ZXt e \ti arcd ti + ti a tcdti )) 

+2c 9 R ef (R r r R t etf + R\Rr f - 2R tr R retf + R bedf (R\ rd + R\ td )) 

+8cio-R 3 + 8cii(2R 2 R tr tr + RR abcd Rabcd) 

+2ci 2 (R 2 (R r r + R t t) + 2RR ab R ab ) + ci 3 (AR abtr (R abtd R d r — R abrd R d t ) + AR trcd R trcg R d 9 

+ (RabcrRef Cr + RabctRef Ct )R abef ) + %Ci A (R r erf R geh fR 9tH t — R^* R geh f R 9 r k t)}, (34) 
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where the q's are given in Appendix. Now integrating over the horizon and dividing by 
V n -i, the entropy density reduces to 

(n-l)l 6 y3 , n _ x 

Combining all of these expressions, the entropy density for quartic quasitopological gravity 
becomes 

rT l A (n-1)/ 2 , 2 „ (n-1)/ 4 A1 , (n-1)/ 6 

4 \ (n — 3) r+ (n — 5 J r_j_ (n — 7) r+ 

A simple method of finding the energy per unit volume V^_i is through the use of first law 



of thermodynamics, dAi = Tds, which gives 18] 



r + /<9s\, (n-l)r?/ , I 2 A )2 / 4 A lS / 6 A l4 / 6 

T ^ — " r + = 77 ^0 + k ~ + A*2« — + /J> 3 k -g- + /i 4 fc 

\ar+ / lo7r \ r+ r+ r° 

The energy density can be written in term of the geometrical mass by use of Eq. (j!5j) as 

M = m. 

16n 

We pause to comment that for k — — 1 it is possible to have negative mass and/or entropy 
for certain values of the couplings. The phenomenon was originally noted for the Einstein 



anti de Sitter number of years ago 



15] • In this case the mass can be negative up to 



a certain extremal value; the entropy is always positive. This situation also can happen in 
Einstein-Gauss-Bonnet gravity. However the quasitopological terms allow for both quantities 
to be negative for certain values of the parameters (a situation that can also occur for third- 



order Lovelock gravity 20]). If the entropy is negative it is not clear what solution is suitable 
as a reference state vacuum solution. 

One approach to treating this problem is to add an overall constant to the entropy such 
that s > {21]]. This approach assumes that Ai > 0, from which a minimal value for r + 
is then obtained using fl36|) . If fa and /14 both vanish, such a minimal value is assured, 
and a minimal value of the entropy is obtained. While having the peculiar feature that the 
entropy vanishes despite the nonvanishing surface area of the minimal r + = r +m j n black 
hole, an ambiguity in the Noether charge approach makes such an assignment possible 
However if 7^ and/or (14 7^ there will be several possible minima for r+, and it is 
no longer completely clear what assignment should be made in order to ensure the entropy 
remain positive. This situation will also hold in all Lovelock theories third-order and higher, 
and we leave a complete treatment of this subject for future study. 
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2i|. 



A general thermodynamic treatment of these black holes can be carried out along lines 
similar to that recently carried out for Lovelock black holes 22]. 



VII. STABILITY OF THE SOLUTIONS 

An investigation of the full stability of the solutions we have obtained is beyond the scope 
of this paper. As a first step, we here consider the stability of the solutions against a class 
of small nonspherical perturbations. 

The metric of a slowly rotating solution in five dimensions may be written as 

dv 2 

ds 2 = -f(r)dt 2 + — — + 2ag(r)h(6)dtd(p + r 2 dtt 2 , (37) 
f( r ) 

where dfl 2 is the metric of a 3-sphere and a is the rotation parameter, which is assumed to 
be small. The first three terms of the action are stable against a nonspherical perturbation, 
while the fourth one (the cubic quasitoplogical term) is stable against small nonspherical 
perturbations [4]. 

Here, we consider the stability of the solutions for the quartic quasitopological term 
against the above class of small nonspherical perturbations. Consider the first scalar term 
in the Lagrangian We find for the slowly rotating metric (1371) that 

R abcd R abe fR efgh R^ h = f 4 + 4 + ^[1 - f(r)] 2 + 0(a 2 ), 

where prime denotes the derivative with respect to the coordinate r. This shows that a 
Lagrangian containing only the first term has no stable solution against a nonspherical 
small perturbation. Upon further investigation we find this happens for all the other terms 
of the quartic quasitopological term given in (jSJ). 

Consider next the Lagrangian (jH]) with the coefficients given in Eq. (J7J). To first order in 
a we find that it reduces to 



*4 = ^ {(1 - f) 3 r 2 f" " 3 [(1 - /) + rff) + Oi 



a 2 ) 



and so the perturbed quartic Lagrangian contains most second order derivatives. Therefore 
the field equation is a second order differential equation to linear order in a, and so the 
spherical solutions are stable. 
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VIII. HOLOGRAPHIC HYDRODYNAMICS 



As a first step in understanding the role of our theory in the context of the AdS/CFT 
correspondence, we compute the ratio of shear viscosity to entropy, rj/s, leaving other sub- 
jects such as the holographic trace anomaly and holographic computation of energy fluxes 
for future study 



The first computations of rj/s from an AdS/CFT perspective appeared in 23| for Ein- 
stein gravity, and leadin g co rrections for strongly coupled N = 4 super- Yang- Mills theory 
subsequently followed (JO). These computations have been carried out for second [ij and 
third-order Lovelock theories js] and quasitopological gravity 0. Further investigations also 



provided increasingly efficient techniques for these calculations [25 
pole method 



26|. Here, we use the 



26|, for the planar class of metrics 



ds 2 = — (^p(r)dt 2 + dx\ + dx\ + dx^) 



l 2 dr 2 



where ip(r) is the root of Eq. ( fi~3l) with /to = 1. 

Employing the transformation 2 = 1 — r^m 1 / 2 , the metric (|38|) becomes 



(38) 



ds' 



(?p(z)dt 2 + dx\ + dx\ + dx^) 



(39) 



l 2 (l-z) vrv 7 1 2 61 4(l-z)V»(z)' 

where ^(2) has a simple zero at the horizon located at 2 = 0. Thus if)(z) may be expanded 

as 

^(z) = $\ + 4 2) z 2 + 4 3) z 3 + 4 4) z 4 + (40) 

where ip^ is the zth derivative of ip(z) at z = 0. Using Eq. (ITS]) and the Taylor expansion 
Hi]) , the expansion coefficients can be obtained as 



,(i) 



to' = -2, 4 2) = 2(1 - 4/i 2 ), = 24(/i 2 - 4/i 2 . + 2/x 3 ), 

^ 4) = 24 [/2 2 - 24/t 2 . + 80/2 2 (A 3 - £2) - 12£ 3 - 16£ 4 ] • 

We perturb the metric (1351 by the shift 

ctej — y dxi + ee~ lwt dxj, 

and we calculate the Lagrangian density. Because of the off-shell perturbatio n ([42]) 
exists a pole at z = in the (otherwise) on-shell action. The shear viscosity is [261 ] 



(41) 

(42) 
there 



V 



-8ttT lim 
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where Res z= oL denotes the residue of the pole in the Lagrangian density, and T is the 
Hawking temperature give in Eq. (USD as T = r+^l 2 )- 1 = m 1 / 4 ^ 2 ) -1 . It is a matter of 
calculation to show that the shear viscosity reduces to 

V = y^{l - 4/z 2 - 36/2 3 (9 - 64// 2 + 128/^ _ 48/} 3 ) 
96 

-—^4 (1491 - 10800/t 2 + 28864/t^ - 6240/t 3 + 10752/} 2 /i3 - 25088/^)}. 

Now, using the fact that the entropy density of the black brane is s = m 3//4 (4Z 3 ) -1 , the ratio 
of shear viscosity to entropy is 

- = ^{l-4/2 2 -36/2 3 (9-64/2 2 + 128/2 2 .-48/2 3 ) 

S 47T 

96 

-— /2 4 (1491 - 10800/i 2 + 28864/}^ - 6240/x 3 + 10752/x 2 /x 3 - 25088/2J) . 

The last term is the effect of the quartic quasitopological term on 77/s. Clearly it can be 
either positive or negative; the investigation of the allowed values of this term will be given 
elsewhere. 



IX. CONCLUDING REMARKS 



We have explicitly constructed the Lagrangian for quartic quasitopological gravity (up 
to a term proportional to the Lovelock term) for all dimensions D > 5 except for D — 8, 
and shown specifically what its black hole solutions are. This is the highest-degree case for 
which it is possible to find explicit solutions. 

It is possible to make some general remarks about quasitopological gravity even though 
the specific Lagrangian has not been found for an arbitrary power K of the curvature. Since 
all derivative terms higher than 2 must be eliminated from the Lagrangian, which itself 
must be linear in the lapse function, it is reasonable to conjecture that the action in the 
spherically symmetric case will be reduce to 

I G = j dtdrN{r) (^ n Y^h^ (43) 

up to terms proportional to the transverse volume V n -\, for K-th order quasitopological 
gravity in (n + 1) dimensions, where ip = l 2 r~~ 2 (k — /) and the fik parameters are rescaled 
coefficients of the fc-th powered curvature term. For a given K this action should be valid 
for all dimensionalities (n+1) larger than 4, except for particular choices where n = 2K — 1. 
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HP 



A similar conjecture was formalized by considering the invariant [12J, |27J 

\f(K) _ rMi^i-WK ( p QijSi p "iffe _ /hi aift ^ ok/^a'A (aa\ 

where ^o 1 '"^ « K is the generalized Kronecker-delta tensor and C„„ a/3 . When rewritten in 
terms of Riemann invariants, factorizes with a common factor of (n — 2fT + 2); for 

n < 2.K" — 2 it vanishes. Taking the action to be of the form 

N { n K) 

N^+Y a {K) W {K 

oK/„ 01/ i or ~ J 3 



2 K (n-2K + 2)" 

3=1 



(45) 



where {W[ , . . . , is a se ^ of linearly independent K-tYi order Weyl invariants, it is 

possible to prove a generalized Birkhoff 's theorem provided olq is an appropriately chosen 
linear combination of the otj coefficients 



All the contractions of K Weyl tensors are 



proportional on spherical/planar/hyperbolic symmetric spacetimes, with no static assump- 
tion required. 

For n > 2K — 1 the invariant J\[^ K > can be expressed as a linear combination of the 
2i^-dimensional Euler density and all conformal invariants. In general this action will yield 
field equations greater than second order, but for the aforementioned linear combination of 
ihe ai coefficients they reduce to second order in the spherical/planar/hyperbolic cases 



121 ] . For n < 2K — 2 the action ( 1431) yields a set of fourth order field equations for an arbi- 



trary metric but a set of vanishing field equations on spherical/planar /hyperbolic symmetric 



spacetimes 



12| 



In view of our results for the quartic case, we propose that there exists i^-th or- 
der quasitopological gravity in any dimension except for n = 2K — 1. For the spheri- 
cal/planar/hyperbolic ansatz ([3]) we conjecture that the nonvanishing action is given in 
131), which yields the field equations 




X>^M =° ( 46 ) 



for the metric function /(r). In general the field equations will be of 4th order, since the 
variation of the action will produce terms proportional to 2nd derivatives of variations of 
metric functions multiplied by powers of the Riemann curvature. Upon integration by parts 
the largest number of derivatives that could act on any term will be 4. 
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This equation has the same form as the corresponding situation in Lovelock gravity [22] , 
the difference being that K < [|] in the Lovelock case, whereas K is not restricted in the 
quasitopological case. The solutions to ( 1461) are given by the solutions to the equation 

X>/ = ^ (47) 

k=0 

which for K > 5 cannot be written explicitly in general. The analysis of the black hole 



solutions for this case completely parallels that of the Lovelock case 22j and we shall not 
repeat it here. 

While the quartic Lagrangian §5§ we have constructed is unique (up to a term proportional 
to the Euler density) insofar as it yields second order differential equations for spherically 
symmetric metrics, its geometrical origins remain somewhat obscure. Since all spherically 
symmetric metrics reduce to effective theories of gravitation in two space-time dimensions, it 
may be that some kind of theorem of principle will single out the choice ([5]) with coefficients 
given in the Appendix. This remains an interesting topic for future study. 

Quasitopological gravity provides a much broader range of parameter space for holo- 
graphic duality. It would be interesting to see what constraints are placed on the entropy- 
to-viscosity ratio for this class of theories, and how they modify their condensed matter 
duals in asymptotically Lifshitz gravity. 
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X. APPENDIX 

Here we present the coefficients of quartic-curvature terms in Eq. ([6]) in n + 1 dimensions. 
Using spherically symmetric metric ([3]), one can show that the Lagrangian ([6]) with the 
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following Cj's 

ci = - (n - 1) (n 7 - 3 n 6 - 29 + 170 n 4 - 349 n 3 + 348 n 2 - 180 n + 36) , 
c 2 = -4 (n - 3) (2 n 6 - 20 n 5 + 65 n 4 - 81 n 3 + 13 n 2 + 45 n - 18) , 
c 3 = -64 (n-1) (3n 2 -8n + 3) (n 2 -3n + 3), 

c 4 = -(n 8 -6n 7 + 12n 6 -22n 5 + 114n 4 - 345n 3 + 468n 2 - 270n + 54), 

c 5 = 16 (n - 1) (10 n 4 - 51 n 3 + 93 n 2 - 72 n + 18) , 

c 6 = - - 32 (n - l) 2 (n - 3) 2 (3 n 2 - 8 n + 3) , 

c 7 = 64 (n-2)(n-l) 2 (4n 3 -18n 2 + 27n-9), 

c 8 = -96 (n - 1) (n -2) (2n 4 - 7n 3 + 4n 2 + 6n - 3) , 

c 9 = 16(n-l) 3 (2n 4 -26n 3 + 93n 2 -117n + 36), 
c 10 = - 31 n 4 + 168 ra 3 - 360 n 2 + 330 n - 90, 
cn = 2(6n 6 -67n 5 + 311n 4 - 742n 3 + 936 72 2 - 576n+126), 
c 12 = 8(7n 5 -47n 4 + 121n 3 -141n 2 + 63n-9), 
ci 3 = 16 n (n - 1) (n - 2) (n - 3) (3 n 2 — 8n + 3) , 

Cl4 = 8 (n - 1) (n 7 - 4n 6 - 15n 5 + 122n 4 - 287 n 3 + 297 n 2 - 126n + 18) , 
reduces to the Lagrangian given in Eq. (II ip. 
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